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Abstract 

Several questions of approximation theory are discussed: 1) can one approxi- 
mate stably in L°° norm /' given approximation fg, \\ fg — f \\l°°< of & n unknown 
smooth function /(x), such that || f'{x) ||l°°< Tn\l 

2) can one approximate an arbitrary / e L 2 (D),D C M n ,ra > 3, is a bounded 
domain, by linear combinations of the products u±U2, where u m S N(L m ),m = 1, 2, 
L m is a formal linear partial differential operator and N(L m ) is the null-space of 
L m in D, N(L m ) := {w : L m w = in D}? 

3) can one approximate an arbitrary L 2 (D) function by an entire function of 
exponential type whose Fourier transform has support in an arbitrary small open 
set? Is there an analytic formula for such an approximation? 

1 Introduction 

In this essay I describe several problems of approximation theory which I have studied 
and which are of interest both because of their mathematical significance and because of 
their importance in applications. 

1.1 

The first question I have posed around 1966. The question is: suppose that f(x) is 
a smooth function, say / G C°°(R), which is T-periodic (just to avoid a discussion of 
its behavior near the boundary of an interval), and which is not known; assume that 
its (^-approximation f s e is known, \\fg — fW^ < 5, where || • is the L°°(R) 

norm. Assume also that ||/'||oo < m i < oo- Can one approximate stably in L°°(R) the 
derivative /', given the above data {5, fg, mi}? 

*key words: stable differentiation, approximation theory, property C, elliptic equations, Runge-type 
theorems, scattering solutions 

tMath subject classification: 41-XX, 30D20, 35R25, 35J10, 35J05, 65M30 



By a possibility of a stable approximation (estimation) I mean the existence of an 
operator Lg, linear or nonlinear, such that 



sup ||W«-/'lloo<f7(<y)->0as<y->0, (1.1) 

/ec°°(R) 
||/-/j||oo<5,||/'||oo<mi 

where T](S) > is some continuous function, 77(0) = 0. Without loss of generality one 
may assume that rj(8) is monotonically growing. 

In 1962-1966 there was growing interest to ill-posed problems. Variational regulariza- 
tion was introduced by D.L. Phillips [2] in 1962 and a year later by A.N. Tikhonov 
in 1963. It was applied in pQ in 1966 to the problem of stable numerical differentiation. 
The method for stable differentiation proposed in PQ was complicated. 

I then proposed and published in 1968 |3] the idea to use a divided difference for 
stable differentiation and to use the stepsize h = h(5) gularization parameter. If 

||/" || < ma, then h(5) = and if one defines (0): 

T f fs(x + h(S)) - f 5 (x - h(S)) [W 

uh '- = 2^ ' m = yw { ] 

then 



\\Lgfg - /'Or)||oc < := e{5). (1.3) 

It turns out that the choice of Lg, made in [Hj, that is, Lg defined in (1.2), is the best 
possible among all linear and nonlinear operators T which approximate f'(x) given the 
information {5, m 2 ,fg}. Namely, if JC(5, rrij) := {/ : / G C J (R),rrij < oo, 1 1 / — /<5 1 1 oo < $}, 
and rrij = \\f^\\oo, then 



inf sup \\Tf s -f\\>y^S. (1.4) 

T feK(5,m 2 ) 

One can find a proof of this result and more general ones in [Oj and various 

applications of these results in [lj - 

The idea of using the stepsize h as a regularization parameter became quite popular 
after the publication of [3 and was used by many authors later. 

In [57] formulas are given for a simultaneous approximation of / and /'. 



1.2 

The second question, that I will discuss, is the following one: can one approximate, with 
an arbitrary accuracy, an arbitrary function f(x) G L 2 (D), or in L P (D) with p > 1, 
by a linear combination of the products u±U2, where u m G N(L m ), m = 1,2, L m is 
a formal linear partial differential operator, and N(L m ) is the null-space of L m in D, 
N(L m ) := {w : L m w = in D}? 
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This question has led me to the notion of property C for a pair of linear formal partial 
differential operators {Li,L 2 }. 

Let us introduce some notations. Let D C R n ,n > 3, be a bounded domain, 
L m u(x) := J2\j\<j m cijm(x)D^u(x), m = 1,2, j is a multiindex, J m > 1 is an integer, 
a,j m (x) are some functions whose smoothness properties we do not specify at the moment, 
DJu = o^lrir- ^ =h + --- + Jn- Define N m := N D (L m ) := {w : L m w = in D}, 
where the equation is understood in the sense of distribution theory. Consider the set of 
products {wiw 2 }, where w m G N m and we use all the products which are well-defined. If 
L m are elliptic operators and a jm (x) G C 7 (IR ?1 ), then by elliptic regularity the functions 
w m G C 7+Jm and therefore the products Wiw 2 are well defined. 

Definition 1.1. A pair {L 1} L 2 } has property C if and only if the set {wiW2\\/ Wm &N m is 
total in L P (D) for some p > 1. 

In other words, if f G L P (D), then 

jjT f(x) Wl w 2 dx = 0, V^ m G AT m | / = 0, (1.5) 

where Vw m G N m means for all w m for which the products W1W2 are well defined. 

Definition 1.2. If the pair {L, L} has property C then we say that the operator L has 
this property. 

From the point of view of approximation theory property C means that any function 
/ G L P (D) can be approximated arbitrarily well in L P (D) norm by a linear combination 
of the set of products W1W2 of the elements of the null-spaces N m . 

For example, if L = V 2 then N( V 2 ) is the set of harmonic functions, and the Laplacian 
has property C if the set of products hih 2 of harmonic functions is total (complete) in 
L p (D). 

The notion of property C has been introduced in ^U]. It was developed and widely 
used in ^U] - [21] ■ It proved to be a very powerful tool for a study of inverse problems 

m - chi, m - eh. 

Using property C the author has proved in 1987 the uniqueness theorem for 3D inverse 
scattering problem with fixed-energy data ^2], [T3*] . |17j . uniqueness theorems for inverse 
problems of geophysics [12], C^!; [IS], and for many other inverse problems [IB]. The 
above problems have been open for several decades. 

1.3 

The third question that I will discuss, deals with approximation by entire functions 
of exponential type. This question is quite simple but the answer was not clear to 
engineers in the fifties, it helped to understand the problem of resolution ability of linear 
instruments [22] , [22] , and later it turned to be useful in tomography [2H] • This question 
in applications is known as spectral extrapolation. 
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To formulate it, let us assume that D C R™ is a known bounded domain, 

7(0:= [ f(x)e*-*dx:=Ff, f(x) E L 2 (D), (1.6) 

and assume that /(£) is known for £ E D, where D is a domain in R^ - . The question is: 
can one recover f(x) from the knowledge of /(£) in D? 

Uniqueness of f(x) with the data {/(£)>£ £ ^} is immediate: /(£) is an entire 
function of exponential type and if /(£) = in D, then, by the analytic continuation, 
/(£) = 0, and therefore f(x) = 0. Is it possible to derive an analytic formula for the 
recovery of f(x) from {f(£), £ E D}1 It turns out that the answer is yes - [2S])- Thus 
we give an analytic formula for inversion of the Fourier transform /(£) of a compactly 
supported function f(x) from a compact set D. 

From the point of view of approximation theory this problem is closely related to the 
problem of approximation of a given function h{£) by entire functions of exponential type 
whose Fourier transform has support inside a given convex region. This region is fixed 
but can be arbitrarily small. 

In sections 2,3 and 4 the above three questions of approximation theory are discussed 
in more detail, some of the results are formulated and some of them are proved. 



2 Stable approximation of the derivative from noisy 
data. 

In this section we formulate an answer to question 1.1. Denote ||/ ( - 1+a ' ) || := mi+ai where 
< a < 1, and 

!l/ (1+a) !lHmU+su P l^-yi . (2.i) 



\x - y\ 

Theorem 2.1. There does not exist an operator T such that 

sup HT/a - /'I,*, < 7/(5)^0 as 5^0, (2.2) 

if 3 = or 3 = 1- There exists such an operator if j > 1. For example, one can take 
T = Lgj where 

r f fs(x + hj(6)) - f s (x - hj(6)) ( 5 y 

and then 

sup \\L s ,jfs - / ||oo < Cj5 j , 1 < j < 2, (2.4) 
j i 

Cj := — r m] . (2.5) 
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Proof. 1. Nonexistence of T for j — and j = 1. 

Let = 0, := - ma(a ~ 2fe) , < x < 2/i. Extend ^(x) on R so that H/^IU = 

suPo<x<2h WfPWJ = 0, 1,2, and set / 2 (x) := -/x(x). Denote (T/i)(0) := 6. 
One has 

||T/ 4 - /{|| > |(T/,)(0) - f[(0)\ = \b- mh\, (2.6) 

and 

\\Tf 5 -f 2 \\ > \b + mh\. (2.7) 
Thus, for j — and j — 1, one has: 

7j := inf sup ||T/,5 — /'|| > inf max [|6 — mh\, \b + mh\] — mh. (2.8) 
Since \\f s - fs\\oc < s = 1, 2, and / 5 = 0, one gets 

ll/.IU = -y < S, 8 = 1,2. (2.9) 

Take 

/i=t/^. (2.10) 



Then m = \\f s \\oo = S,m 1 = U/J^, = V25m, so (2.8) yields 

7o = \ / 25m — > oo as m — > oo, (2-H) 



and (2.2) does not hold if j = 0. 
If j = 1, then (2.8) yields 



7i 



= V2d^m = mi >0, (2.12) 



and, again, (2.2) does not hold if j = 1. 
2. Existence of T for j > 1. 

If j > 1, then the operator T = Lsj, defined in (2.3) yields estimate (2.4), so (2.2) 

j-i 

holds with t)(5) = Cj5 i and Cj is defined in (2.5). Indeed, 

\\Lsjfs - /'Hoc < \\LsMs ~ /)l|oc + \\L SJ f - /'Hoc < I + m^'" 1 . (2.13) 

Minimizing the right-hand side of (2.13) with respect to h > for a fixed 5 > 0, one gets 
(2.4) and (2.5). 

Theorem 2.1 is proved. □ 

Remark 2.1. If j = 2, then m 2 = m, where m is the number introduced in the beginning 
of the proof of Theorem 2.1, and using the Taylor formula one can get a better estimate 
in the right-hand side of (2.13) for j = 2, namely \\L$ >2 fs — /'||t» < f + ^^f^. Minimizing 

with respect to h, one gets h(5) = \J~^ and min^o (f + In 2 1 ) : = s(S) = y / 25m 2 . Now 
(2.8), with m = m 2 , yields 

7 2 > y/25m 2 = £(5). (2.14) 
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Thus, we have obtained: 
Corollary 2.1. Among all linear and nonlinear operators T , the operator Tf = Lgf := 

"<<» = sM- 



^ , ' r ""'' 1 - "' " >s ^\ h{5) = J—, yields the best approximation of f, f E fC(S, m 2 ), and 



72 :=inf sup \\Tfs-f'\\ 00 = e(5):= y f 26n7 2 . (2.15) 

T feK(6,m 2 ) 

Proof. We have proved that 72 > s(S). If T = L$ then \\Lsfs — f'\\<x < £ (.$)i as follows 
from the Taylor's formula: if h = \ —, then 

l|W.-/'ll<^ + ! ^ = ^), (2.16) 
so 72 = e(5). □ 



3 Property C 
3.1 

In the introduction we have defined property C for PDE. Is this property generic or is it 
an exceptional one? 

Let us show that this property is generic: a linear formal partial differential operator 
with constant coefficients, in general, has property C.In particular, the operators L = 
V 2 , L = id t — V 2 , L = d t — V 2 , and L = <9 2 — V 2 , all have property C. 

A necessary and sufficient condition for a pair {L 1; L 2 } of partial differential operators 
to have property C was found in [TU] and |2H1 (see also [TS]). 

Let us formulate this condition and use it to check that the four operators, mentioned 
above, have property C. 

Let L m u := <\<j m a jm,D^u(x), m = l,2,a jm = const, x G M n , n > 2, J m > 1. 
Define the algebraic varieties 

£ m := {z : z e C n , L m (z) := ^ a jm z J =0}, m = 1, 2. 

b'|<Jm 

Definition 3.1. Let us write L\ \ C 2 if and only if there exist at least one point E C\ 
and at least one point z^ e £ 2 , such that the tangent spaces in C n T m to C m at the points 
y m ', m — 1,2, are transversal, that is T\\T 2 . 

Remark 3.2. A pair {£\, C 2 } fails to have the property C\ \ C 2 if and only if C\ U C 2 
is a union of parallel hyperplanes in C n . 

Theorem 3.3. A pair {Li,L 2 } of formal linear partial differential operators with con- 
stant coefficients has property C if and only if Ci $ C 2 . 
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Example 3.4 Let L = V 2 , then C = {z : z G C n , z\ + z\ + . . . z 2 n = 0}. It is clear 
that the tangent spaces to C at the points (1,0,0) and (0, 1,0) are transversal. So the 
Laplacian L = V 2 does have property C (see Definition 1.2 in section 1). In other words, 
given an arbitrary bounded domain D C M n and an arbitrary function f(x) G L P (D), 
p > 1, for example, p = 2, one can approximate /(x) in the norm of L P (D) by linear 
combinations of the product h x h 2 of harmonic in L 2 (D) functions. 

Example 3.5 One can check similarly that the operators d t — V 2 , id t — V 2 and d 2 — V 2 
have property C . 

Proof of Theorem 3.3 We prove only the sufficiency and refer to JH] for the necessity. 
Assume that d |f C 2 . Note that e x ' z G N(L m ) := N m if and only if L m (z) = 0, z G C n , 
that is z G C m . Suppose 



The function F(z), defined in (3.1), is entire. It vanishes identically if it vanishes on an 
open set in C" (or IR n ). The set {zi + z 2 }^z m eN m contains a ball in C n if (and only if) 
Ci Indeed, if Z\ runs though the set L\ PI B(z^',r), where B{z { - l \r) := {z : z G 

C n , l^^ 1 )— 2 1 < r}, and z 2 runs through the set C2^B{z^ 2 \ r), then, for all sufficiently small 
r > 0, the set {z\ + 22} contains a small ball p), where C = z^ 1 ^ + z^ 2 \ and p > is a 
sufficiently small number. To see this, note that Ti has a basis hi, . . . h n _i, which contains 
n — 1 linearly independent vectors of C n , and T 2 has a basis such that at least one of its 
vectors, call it h n , has a non-zero projection onto the normal to T 1; so that {hi, . . . , h n } 
are n linearly independent vectors in C n . Their linear combinations fill in a ball B{C,p). 
Since the vectors in T m approximate well the vectors in C m fl B(z^ m \r) if r > is 
sufficiently small, the set {z 1 + z 2 }\/ Zme c m contains a ball B((,p) if p > is sufficiently 
small. Therefore, condition (3.1) implies F(z) = 0, and so f(x) = 0. This means that the 
pair {Li, L 2 } has property C. We have proved: {Ci || C 2 } ^> {{-^1, L 2 } has property C}. 
□ 

How does one prove property C for a pair {Li, L 2 } := {V 2 + A; 2 — qi(x), V 2 + k 2 — q 2 (x)} 
of the Schrodinger operators, where k = const > and q m (x) G L 2 oc (M. n ) are some real- 
valued, compactly supported functions? 

One way to do it [18J is to use the existence of the elements ip m G N m := {w : 
[V 2 + k 2 — q m (x)]w = in IR n } which are of the form 




then 




(3.1) 



i) m {x, 0) = e ike - x [l + R m {x, 9)}, k > 0, (3.2) 

where 9 G M := {9 : 9 G C n , 9 ■ 9 = 1}. Here 9 ■ w := YT j= i ®i w h ( note th &t there is no 
complex conjugation above Wj), the variety M is noncompact, and [TH] 
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00, m = 1, 2, 



(3.3) 
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where c = const > does not depend on 9, c depends on D and on ||9||i°°(s a ), where 
q = q, q = for \x\ > a, and D C M. n is an arbitrary bounded domain. Also 

\\Rm(x,e)\\ L 2 (D) <-^- 0eM, |0|->oo, m = l,2. (3.4) 

It is easy to check that for any £ G M n , n > 3, and any > 0, one can find (many) 6\ 
and 6*2 such that 

k(6 1 + 6 2 )=Z, — »• oo, 6 u 9 2 eM, n>3. (3.5) 

Therefore, using (3.5) and (3.3), one gets: 

lim fafo = e ifx , (3.6) 

|9i|-»oo 

Since the set {e^ ,a: }v?eR n is total in L P (D), it follows that the pair {Li,L 2 } of the 
Schrodinger operators under the above assumptions does have property C. 



3.2 

Consider the following problem of approximation theory 

Let k = 1 (without loss of generality), a G S 2 (the unit sphere in IR 3 ), and u := u(x, a) 
be the scattering solution that is, an element of N(V 2 + l—q(x)) which solves the problem: 

[V 2 + 1 - q(x)] u = in IR 3 , (3.5) 

u = exp(ia ■ x) + A(a', a)—— + o [ — ) , \x\ — > oo, a \= - — r , (3.6) 

\x\ \\ X \J \ x \ 

where a G S 2 is given. 

Let w G A^(V 2 + 1 — g(x)) := N(L) be arbitrary, G i^; 2 c , ^ is the Sobolev space. 
The problem is: is it possible to approximate w in L 2 (D) with an arbitrary accuracy by a 
linear combination of the scattering solutions u{x, a)l In other words, given an arbitrary 
small number e > and an arbitrary fixed, bounded, homeomorphic to a ball, Lipschitz 
domain D C M n , can one find v e (a) G L 2 (S 2 ) such that 

||it>— / u(x, a)v £ (a)da\\ L 2( D -) < el (3.7) 
Js 2 

If yes, what is the behavior of ||t> e (aO ||l 2 (s2) as e — > 0, if w = vf)(x,9) where ip is the 
special solution (3.2) - (3.3), 6 G M, Im9 ^ 0? 

The answer to the first question is yes. A proof [TH] can go as follows. If (3.7) is false, 
then one may assume that w G N(L) is such that 

w(^J u(x, a)v(a)daj dx = Wv G L 2 (S 2 ). (3.8) 
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This implies 

/ Wu(x, a)dx = Va G S 2 . 
Jd 

From (3.9) and formula (5) on p. 46 in ^H] one concludes: 

v( y ) ■= / wG{x,y)dx = G D' :— M?\D, 
Jd 

where G(x, y) is the Green function of L: 

[V 2 + 1 - q(x)] G(x, y) = -8(x, y) in M 3 , 
dG 



lim 

™ >\x\=r 



iG 



ds = 0. 



d\x\ 

Note, that in [|THJ , P- 46] formula (5) is: 

e i\y\ 

G(x,y,k) = - — —u(x,a 
4tt I ?/ 1 

From (3.10) and (3.11) it follows that 

[V 2 + 1 - q(x)] v{x) = -w(x) in D, 
v — v n = on S := dD, 







•a- 


1 2/1 





1 1/1 



-a. 



(3.9) 

(3.10) 

(3.11) 
(3.12) 

(3.13) 

(3.14) 
(3.15) 



loc 



where iV is the outer unit normal to S, and (3.15) holds because v = in D' and v G H } 
by elliptic regularity if q G Lf oc . Multiply (3.14) by w, integrate over D and then, on the 
left, by parts, using (3.15), and get: 



M dx = 0. 



(3.16) 



D 



Thus, w = and (3.7) is proved. □ 
Let us prove that ||u e (a)||x,2( S 2) — ► oo as e — > 0. Assuming ||i> e (a)|| < c We G (0,£q), 
where So > is some number, one can select a weakly convergent in L 2 (S 2 ) sequence 
v n (a) — > w(a), n — >■ oo, and pass to the limit in (3.7) with w = ip(x, 9), to get 



U{x,9) 



S 2 



u(x,a)v(cx)da\\L2(D) = 0- 



(3.17) 



The function U(x) := J 52 u{x, a)v(a)da G A^(L) and ||?7(x)||x,oo(ii») < oo (since sup xgR n Qg5 2 
\u(x, a)\ < oo). Since (3.17) implies that ip(x, 9) = U(x) in D, and both ijj{x, 9) and U(x) 
solve equation (3.5), the unique continuation principle for the solutions of the elliptic 
equation (3.5) implies U(x) = tp(x,9) in IR 3 . This is a contradiction since ip(x,9) grows 
exponentially as |x| — >• oo in certain directions because lm9 ^ (see formula (3.2)). 
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We have proved that if w = if)(x,9), ImQ ^ 0, then ||^ e (a) ||l 2 (s 2 ) c» as e 0, 
where v e {a) is the function from (3.7). 

For example, if q(x) = and k = 1 then i/j(x, Q) = e ld x and u(x, a) = e ia ' x . So, if 
ImQ ^ 0,9 e M, and \\e %6 ' x — J s2 v £ (a)e ta ' ||l 2 (d) < £, then ||z/ e (a)|| L 2( 5 2) — > oo as e — > 0. 
It is interesting to estimate the rate of growth of 



This is done in £U (and |22|). 

Property C for ordinary differential equations is defined, proved and applied to many 
inverse problems in [T§] . 

4 Approximation by entire functions of exponential 



Let f(x) G L 2 (B a ), B a := {x;x G 1R™, |x| < a}, a > is a fixed number, f(x) = for 
|x| > a, and 



Assume that /(^) is known for all £ G D C R^, where D is a (bounded) domain. 

The problem is to find /(£) for all £ G R?, (this is called spectral extrapolation), 
or, equivalently, to find f(x) (this is called inversion of the Fourier transform /(£) of a 
compactly supported function f(x), supp/(x) C B a , from a compact D). 

In applications the above problem is also of interest in the case when D is not nec- 
essarily bounded. For example, in tomography D may be a union of two infinite cones 
(the limited-angle data). 

In the fifties and sixties there was an extensive discussion in the literature concerning 
the resolution ability of linear instruments. According to the theory of the formation of 
optical images, the image of a bright point, which one obtains when the light, issued by 
this point, is diffracted on a circular hole in a plane screen, is the Fourier transform /(£) 
of the function f(x) describing the light distribution on the circular hole B a , the two- 
dimensional ball. This Fourier transform is an entire function of exponential type. One 
says that the resolution ability of a linear instrument (system) can be increased without 
a limit if the Fourier transform of the function f(x), describing the light distribution 
on B a , can approximate the delta-function with an arbitrary accuracy. The above 
definition is not very precise, but it is usual in applications and can be made precise: it 
is sufficient to specify the metric in which the delta-function is approximated. For our 
purposes, let us take a delta- type sequence Sj (£) of continuous functions which is defined 




|Ka)|| L 2 (52) . 



type. 




(4.1) 
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lim / «,.(flde 



by the requirements: | J5 <5j(£)d£| ^ c > where the constant c does not depend on D and 
j, and 

' 1 if e D, 
if £ 5. 

TTie approximation problem is: can one approximate an arbitrary continuous (or 
) function g(£) in an arbitrary fixed bounded domain D C M.^ by an entire function 

of exponential type /(£) = j B f(x)e l ^' x dx, where a > is an arbitrary small number? 

The engineers discussed this question in a different form: can one transmit with an 
arbitrary accuracy a high-frequency signal g(£) by using low-frequency signals /(£)? The 
smallness of a means that the "spectrum" f(x) of the signal /(£) contains only "low 
spatial frequencies" . 

From the mathematical point of view the answer is nearly obvious: yes. The proof is 
very simple: if an approximation with an arbitrary accuracy were impossible, then 

= J_g(Ji) QT e* x f(x)dxj di V/ e L\B a ). 



This implies the relation 



= / g(Oe* x d£ Vx e B a . 

J D 



Since D is a bounded domain, the integral above is an entire function of x G C n which 
vanishes in a ball B a . Therefore this function vanishes identically and consequently 
g(£) = 0. This contradiction proves that the approximation of an arbitrary g(£) £ L 2 (D) 
by the entire functions f(£) = f B f(x)e l ^ x dx is possible with an arbitrary accuracy in 
L 2 (D). 

Now let us turn to another question: how does one derive an analytic formula for 
finding f(x) if /(£) is given in D7 

In other words, how does one invert analytically the Fourier transform /(£) of a 
compactly supported function f(x), supp / C B a , from a compact D? 

We discuss this question below, but first let us discuss the notion of apodization, 
which was a hot topic at the end of the sixties. Apodization is a method to increase the 
resolution ability of a linear optical system (instrument) by putting a suitable mask on the 
outer pupil of the instrument. Mathematically one deals with an approximation problem: 
by choosing a mask g(x), which transforms the function f(x) on the outer pupil of the 
instrument into a function g(x)f(x), one wishes to change the image /(£) on the image 
plane to an image <5j(£) which is close to the delta-function and therefore increase 
the resolution ability of the instrument. That the resolution ability can be increased 
without a limit (only in the absence of noise!) follows from the above argument: one 

can choose g(x) so that g(x)f(x) will approximate arbitrarily accurately 5j(£), which, in 
turn, approximates 6(£) arbitrarily accurately. 
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This conclusion contradicts to the usual intuitive idea according to which one cannot 
resolve details smaller than the wavelength. 

In fact, if there is no noise, one can, in principle, increase resolution ability without 
a limit (superdirectivity in the antenna theory), but since the noise is always present, in 
practice there is a limit to the possible increase of the resolution ability. 

Let us turn to the analytic formula for the approximation by entire functions and for 
the inversion of the Fourier transform of a compactly supported function from a compact 

D - 

Multiply (4.1) by (27r)~ n 8j(£,)e~^' x and integrate over D to get 

m = [ fWit? - V)dy = J^L 7me-^ (4.2) 
JB a y^ 71 ) Jd 

where 5j(£) is the Fourier transform of Sj(x). 

Let us choose 5j(x) so that it will be a delta- type sequence (in the sense defined 
above). In this case fj(x) approximates f(x) arbitrarily accurately: 

\\f-Mx)\\ ^Oasj^oo, (4.3) 

where the norm || • || is L 2 (B a ) norm if / G L 2 (B a ), and C(5 a )-norm if / G C(B a ). 
If ll/llci(Ba) < m U then 

||/-/i(^)||c (Bo )<cm 1 r | . (4.4) 
The conclusions (4.3) and (4.4) hold if, for example, 

6j(x) := Pj(\x\ 2 ) (F^h) (x), (4.5) 



where 

PAr\ ■-( 



P i( r ) := ( 7^2 ) ( 1 - TT2 ) , < r < a, a, > a, (4.6) 



h(OeC™(D), *Jh{t)d£=l. (4.7) 

Theorem 4.1. If (4-5)-(4-7) hold, then the sequence fj(x), defined in (4-2), satisfies 
(4.3) and (44). 

Thus, formula (4.2): 

f(x) = T~ x [J(05(0] » ( 4 - 8 ) 

where 6j(£) := j u 8j(x), and Sj(x) is defined by formulas (4.5)-(4.7), is an inversion formula 
for the Fourier transform /(£) of a compactly supported function f(x) from a compact 
D in the sense (4.3). A proof of a Theorem similar to 4.1 has been originally published 
in |23- 
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In j22], apodization theory and resolution ability are discussed. In j2H] a one- 
dimensional analog of Theorem 4.1 is given. In this analog one can choose analytically 
explicitly a function similar to h(£). 

Proof of Theorem 4.1 is given in [j2S], PP- 260-263]. 

Let us explain a possible application of Theorem 4.1 to the limited-angle data in 
tomography. The problem is: let f(a,p) : = jj f(x)ds, where Z ajP is the Radon transform 

of f(x). Given f(a,p) for all p G R and all a G K, where K is an open proper subset of 
S 2 , find f(x), assuming supp / C B a . It is well known [22], that 

/oo 
f(a,p)e ipt dp = f(ta), teR, tct := £. 
-oo 

Therefore, if one knows f(a,p) for all a G if and all p G R, then one knows /(£) for 
all £ in a cone K X R. Now Theorem 4.1 is applicable for finding /(x) given f(a,p) for 
a G if and pel. 
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